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(3 ■ ABSTRACT 



Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi- 
wavelength information of the CDFS to assess this issue. We select the AGNs from the 3 Ms XMM-Newton survey and measure 
the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 jim, predominantly from deep 
, 100 yL/m and 160 fim Herschel observations, but also from Spitzer MIPS-70 fim. Star-formation rates are obtained from spectral energy 

r \ , distribution fits, identifying and subtracting an AGN component. Our sample consists of sources in the z ~ 0.5 - 4 redshift range, 

^ ■ with star-formation rates SFR ^ 10^ - 10^ Moyr"^ and stellar masses ^ 10^^ - 10^^^ M©. We divide the star-formation rates 

by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between 
the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding 
Lx - 10^^ erg s~^ and redshifts z > I. We do not find evidence for such a correlation for lower luminosity systems or those at lower 
redshifts, consistent with previous studies. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption 
derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather 
than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming 
galaxies at the same redshift (the so-called "main sequence") we find that the AGNs reside mostly in main-sequence and starburst 
hosts, reflecting the AGN - sSFR connection; however the infrared selection might bias this result. Limiting our analysis to the 
highest X-ray luminosity AGNs (X-ray QSOs with Lx > 10"^"^ erg s~^), we find that the highest-redshift QSOs (with z > 2) reside 
predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of 
QSOs at z ~ 1.5 with specific star-formation rates compatible with the main- sequence, or even in the "quiescent" region. Finally, 
we test the reliability of the colour-magnitude diagram (plotting the rest-frame optical colours against the stellar mass) in assessing 
host properties, and find a significant correlation between rest-frame colour (without any correction for AGN contribution or dust 
extinction) and sSFR excess relative to the "main sequence" at a given redshift. This means that the most "starbursty" objects have 
the bluest rest-frame colours. 



Key words. Galaxies: active - Galaxies: Seyfert - Galaxies: statistics - Galaxies: star formation - X-rays: galaxies - Infrared: 
galaxies 



1. Introduction 

One of the most significant observations of modern-day astro- 
physics is the evidence that the mass of the super-massive black 
hole (SMBH) in the centre of any galaxy is correlated to the 
prop erties of its bulge , parametrised by the spheroid luminosity 
(e.g.lMagorrian et al.. 1998), or the spheroid velocity dispersion 
(e.g. Herrarese & Merritt , 2QO 0). Thi s relati on has a small in- 
trinsic dispersion (e.g. Giiltekin et al., 2009) which implies an 
evolutionary connection between the SMBH and the spheroid. 
The mechanisms that build the super-massive black hole and the 
bulge of the galaxy are an active galactic nucleus (AGN) and 
star-formation or possibly merging episodes, respectively. There 
is additional evidence that the space density of AGNs and cos- 
mic star formation ha ve similar redshift evolution, at least up 
to red shifts z ~ 2 (e.g. IChapman et al.L 120051: iMerloni & Heinzl. 
I2008h . 

The coeval growth of the SMBH and the host galaxy 
implies some causal con nection between the AGN and star- 
formation properties (see [Alexander & Hickoxi 1201 2i for a re- 
view). Theoretical and semi-analytical models of galaxy evolu- 
tion through mergers assume s uch a connectio n, wher e AGN 
feedback (e.g. Hopk ins et all 120061: Di MatteoetaD, l2008l) 
plays a catalytic role. After the SMBH has grown sufficiently 
massive, the outflows driven by the radiation pressure of the 
AGN have enough energy to disru pt the cold gas supply which 
sustains the star formation (e.g. Springel et al.[ 120051: iKingl 
I2005h . giving rise to the SMBH-bulge relation. The gas sup- 
ply for both the AGN and the star formation is often thought to 
come from the galaxy mergers, which are ideal mechanisms for 
removing angular momentum from the participant g alaxies and 
funnelling gas to the central kpc region (e.g. iDi Matteo et al.L 
^2005; Barnes & Hernquist, 1996). 

There is, however, growing evidence that a significant part 
of galaxy evolution takes place in secularly evolving systems. 
There is a well-defined relation between the star-formation 
rate and the stellar mass in local star -forming systems (see 
e.g. iBrinchmann et al 1 120041: ISalim et al... i2007) which defines 
the so-called "main sequence" of star formation This rela- 
tion is also found in hig her redshift galaxies (e.g. 'Elbaz et all 
120071: iDaddi et al.l. l2007h with a redshift-dependent normalisa- 
tion. It is also observed that more signs of recent merging ac- 
tivity are found in the morphology of starbursts (defined as 
star-forming galaxies with star-formation rates higher than the 
main sequence) than in normal (main- sequence) star-forming 
galaxies (Kartaltepe et al., 2012). Mapping the star-formation 
in high-redshift (z ~ 1-3) galaxies u sing integral field spec- 
troscopy, Forste r Schreiber et al.l (l2009l) find that about one third 
of those star-forming galaxies have rotation-dominated kinemat- 
ics sh owing no signs of mergers. Moreover, Rodighiero et alj 
(|2()ll|) have shown that ~ 90% of the star-formation den- 
sity at z ~ 1-3 takes place in the main- sequence galaxies. 
The hosts of AGN s do not seem to significantly deviate from 
this main sequence ([Mullaney et al.Ll2012al:ISantini et al.Ll2012l) . 
SimilarlyjGrogin et al.l (120051) found no apparent connection be- 
tween mergers and AGN activi ty at redshifts 0.4 < z ^ 1.3, a re- 
sult which is also confirmed by lCisternas et al.l (1201 lb in a simi- 
lar redshift range (0.3 < z < 1.0). and bv lKocevski et aP (l2012h 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

This work is based on observations obtained with XMM-Newton, 
an ESA science mission with instruments and contributions directly 
funded by ESA Member States and the USA (NASA). 



at higher redshifts (1.5 < z < 2.5). In this case, gravitational 
instabilities of the system may cause the transfer of material 
to the centre through the formati on of bars and pseudo -bulges 
(Kormendv&Kennicutt, 2004 ). iHopkins & Quataeri (l201Qh 
and Diamond- Stan ic & Rieke (i2012h connected the black-hole 
accretion to the nuclear star-formation. In this study, we expand 
the search for an AGN-host connection to higher redshifts. 

Observationally the identification of a connection between 
the star-formation and accretion rates is challenging, especially 
at high redshifts. The most efficient way is to isolate the char- 
acteristic emission bands of both processes, namely the hard 
X-ray emission from the hot corona of the AGN and the far- 
infrared emission from cold dust heated by the UV radiation of 
massive young stars or radio synchrotron emission from elec- 
trons accelerated in supernova explosions. Previous studies us- 
ing those in dicators in deep fiel ds have shown hints of a cor- 
relation (e.g. iTrichas et al.L '2009), which is m ore prominent in 
AGNs with higher luminosities and red shifts (iMullaney et al.L 
l2010l:lLutz e"talll2010l:IShao et al.Ll2010l) . These results argue in 
favour of diff'erent mechanisms, secular evolution and evolution 
through mergers, which take place at lower and higher redshifts 
(or low er and higher luminosities), respectively. Mulla ney et all 
(|2012a|) caution about the eff'ects of both the X-ray (i.e. AGN) 
and the infrared (i.e. star formation) luminosities increasing with 
redshift, which could mimic a correlation between those values, 
especially in samples spanning orders of magnitudes in both Lx 
and LiR, and find no clear signs of a correlation between Lx and 
LiR in moderat e luminosity AGNs (Lx = 10"^^ - 10"^"^ erg s"^). 
More recently, iMuUaney et al.l (l2012bl) do find hints of coeval 
growth of the super-massive black hole a nd the host galax y sug- 
gesting a causal connection (see also Rosa rio et al .l,l201 2| ). 

In this paper we use the deepest observations from XMM- 
Newton and Herschel, combined with Chandra positions and 
deep multi- wavelength data in the CDFS to investigate the AGN- 
host connection, expanding to the less well-sampled region of 
high X-ray luminosities (Lx > 10^^ erg s"^) and redshifts (z > 
2.5). We exploit the multi- wavelength information implement- 
ing an accurate SED decomposition technique to disentangle 
the AGN and star-formation signals in the optical and infrared 
bands, and therefore obtain unbiased star-formation rates for the 
AGN sample. We also make use of accurate XMM-Newton spec- 
tra from the deepest 3 Ms observation for the first time, to inves- 
tigate the nature of the AGN - star-formation relation. 



2. Data 

2.1. X-rays 

Our X-ray data come from the 3 Ms CDFS XM M-Newton sur- 
vey. I nitial results of the survey are presented in lComastri et al.l 
('2011), and details on the data analysis and source detection 
will be presented in Ranalli et al. (in preparation). Briefly, the 
bulk of the X-ray observations were made between July 2008 
and March 2010, and have been combined with archival data 
taken between July 2001 and January 2002, using a single point- 
ing, and covering a total area of 30 x 35 arcmin, centred at the 
Chandra pointing of the CDFS. The total integration time of use- 
ful data is ^ 2.82 Ms. Standard XMM-Newton software and pro- 
cedures were implemented for the analysis of the data, yielding 
a point- spread function (PSF) FWHM of ^ 10.5 arcsec, which 
does not show a significant variation with the off'-axis angle. The 
XMM-CDFS main catalogue contains 337 sources detected in 
the 2-10 keV band with a > 4(T significance, plus a list of 74 
supplementary sources (detected with PWXDetect, but not with 
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EMLDetect), down to a flux limit of ~ 6.6 x 10"^^ erg s"^ cm"^. 
X-ray spectra are produced for 169 sources from both lists, de- 
tected with a significance above Scr and a flux limit of ~ 2x 10"^^ 
erg s"^ cm"^. The spectra have been fitted in XSPEC with a sim- 
ple baseline model of an absorbed power-law and the addition, 
if necessary, of a soft excess component and an Fe Ka line. 



2.2. Optical - near-Infrared 

The area around the CDFS is one of the best observed areas in 
the sky, with a wealth of data. In this work, for the identification 
of our sources in the opt ical and near-I R wavelengths we use 
the MUS YC catalogu es of'Ga wiser et al.1 (12006) and T avlor et al.1 
(l2009h . iGawiser et al . (2006) present the o ptical survey of the 
extended CDFS (see iLehmer et al.l. I2005L hereafter E-CDFS) 
with the MOSAIC II camera of the 4-m CTIO telescope, using 
a BVRIz' filter set. The source extraction is done using a com- 
bined BV R image and the catalogue is complete to Rab = 25. 
iTavlor et al. ( 2001 ) comb ine a large set of optical data, including 
the Gawiser et al 1(12006) data- set, with near-IR data, primarily 
from the ISPI instrument on the CTIO telescope. The catalogue 
contains sources detected in the K band down to a 5 cr limit of 
^AB = 22 and includes photometry in the UUssBVRIz' JHK 
bands. 



2.3. Mid-Infrared 

The entire MUSYC area has been imaged with Spitzer-IRAC in 
four bands, 3.6, 4.5, 5.8, and S.Oyum. The central region is im- 
aged as part of the GOODS survey, and these data are combined 
with more recent observations o f the w ider E-CDFS area in the 
SIMPLE survey (Da men et al .l 1201 ll). The combined data- set 
has a 5 (T magnitude limit of [3.6 yum]AB = 23.86, while the 3 cr 
magnitude limit of the central GOODS region is [3.6yum]AB = 
26.15. 

The GOODS area in the centre of the CDFS has been imaged 
with Spitzer-MlPS in the 24 jim band with a 5 cr flux density 
limit of 30 yuJy. A much wider area, including the entire E-CDFS 
was imaged as part of the FIDEL legacy program (PI: Dickinson; 
description in Magnelli et al., 2009) with a 5 cr flux density limit 
of 70yuJy; we use a combination of the two data-sets for this 
work. 



2.4. Far-Infrared - sub-mm 

The entire E-CDFS region has been imaged with Spitzer MIPS 
in the 70 yum band as part of the FIDEL survey. For the inner 
part of the field we also use the c ombination of obs ervations 
from the GOODS-Herschel survey (lElbaz et al.L l201lh a nd the 
PACS Evolutionary Probes programme (iLutz et al nmB. This 
combination provides the deepest survey of Herschel using the 
PACS instrument (Poglitsch et al, 2010) in both the 100 and 
the 160yum bands, with a total integration time of more than 
400 hours. Because GOODS-Herschel observations cover only a 
13' X 1 r field inside the GOODS area, the combined GOODSH- 
PEP observation has inhomogeneous coverage of the GOODS-S 
field, with the GOODS-S outskirt being 2 times shallower than 
the inner deep area. The data reduction and image construction 
procedures for the FI DEL and the GOODS-H erschel surveys are 
described in detail in Magnelli et al.' ("2009") and Magnelli et al. 
(in prep; but see also Elbaz et al. 201 1 and Lutz et al. 201 1), re- 
spectively. For the source identification and flux density determi- 
nation, all images (MIPS-70 and PACS) were treated in a consis- 



tent way: MIPS and PACS flux densities were derived with a PSF 
fitting analysis, guided using the position of sources detected 
in the deep MIPS-24 observat ions described in Sect. l23l This 
method, presented in detail in lMagnelli et al.l (l2009i 1201 ll) . has 
the advantage that it deals with a large part of the blending issues 
encountered in dense fields and provides a straightforward asso- 
ciation between MIPS and PACS sources. This MIPS-24-guided 
extraction is also very reliable for the purpose of this study, be- 
cause in the GOODS-S field the MIPS-24 observations are deep 
enou gh to contain al l the AGNs of the M I PS-70 and PACS im- 
ages (iMagnelh etall l201ll: iMagdis et all l201lh . The flux den- 
sity limits of the MIPS-70 catalogue used is 2.5 mJy (6cr). For 
the PACS 100 and 160yum, flux density limits of our catalogues 
is 0.6 and 1 .2 mJy (3 cr) in the 1 3' x 1 r inner part of the GOODS- 
S field and 1.2 and 2.4 mJy (3 cr) in the outskirt of the field. AH 
these values include confusion noise. For the sub-mm part of the 
spectrum, we also use the 870 m LA BOCA and 1 .1mm AzTEC 
catalogues of lWeiB et al.1 (12009) and IScottetaD (l2010h . which 
reach depths of 3.5 mJy and 1.4 mJy at the 3.7 and 3.5 cr levels, 
respectively. 



2.5. Radio 

The E-CDFS has been observed with the VLA in two bands (20 
and 6 cm) and the catalogues a re presented in lKellermann et al] 
(l2008h and lMiller et al.1 (l2008h . the former presenting both the 
20 and 6 cm results, and the latter presenting the deeper 20 cm 
catalogue. The 5 cr flux density limit of the survey near the cen- 
tre of the field is 43 yuJy and 55 //Jy, at 20 cm and 6 cm, respec- 
tively. For this work we also chec k the much wider a nd shal- 
lower ATCA 20 cm observations of ^ Norris et al.l (l2006l) , but we 
do not find any new identifications of X-ray sources, however we 
do find some unique spectroscop ic re dshift measurements from 
their follow-u p program; see Sect.l2.6[ We also use the VLBI cat- 
alogue of Mid delberg et al.l (l201ll) to identify any high surface- 
brightness VLBI cores among the radio detections, suggestive of 
high surface-brightness AGN cores. 



2.6. Redshifts 

There are a number of spectroscopic campaigns of the CDFS 
and the E-CDFS available in the literature. For the purposes 
of this paper w e use spectro sc opic redshifts from the fol- 



lowing works: IBalestra et"aD (|2010|) : ICasev et aP 
ICooper et al] (l20n|)dKriek et al.' ('2008^:'Le Fevre et"^ 



2011) 



2004) 



Le Fevre et al.1 d2005h: [M ignoH et al. ( 20 051): iN^ris et al 



( 2006); Ravikumar et al. ( 2007); Silverman et al. (2010) 
Szokolv et al. (2004): Tavlor et al. (2009): Treister et a l. (2009), 
vanderWel et al.1 (l2005b and lVanzeUa et al.l (l2008h . For sources 



which have no spectroscopic redshift determination we use pho- 
tometric redshift estimates from Cardamone et al. (1201 Oah who 
use up to 32 optical and infrared bands, including 18 medium 
narrow-band filters, for ^VT^-detected sources in the E-CDFS. In 
cases where the redshift is not available in the lCardamone et all 
(2010a) catalogue, or it is flagged as low-quality, we use the 
photometric redshifts of Tavlor et al. ( 2009) using 10 bands on 
^-selected sources, iRaff'erty et al.l (2011) who use publically 
available photometric catalogues to determine the ph otometric 
redshifts of E-CDFS sources, and iLuo et al" I (12010) who use 
up to 35 bands from public catalogues to derive redshifts of 
counterparts of Chandra 2-Ms sources. The typical scatter of 
the photometric redshifts is Az/(l + z) ~ 0.01, and using this 
value, we estimate that the induced uncertainty in the infrared 
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Fig. 1: Spatial limits of the different surveys used in this work. 
The grey-scale image is the combined 2-10 keV 3 Ms image of 
the XMM-Newton observations, and the regions are the Herschel 
area (small rectangle), the 4 Ms Chandra area (solid circle), 
the XMM-Newton area used (dashed circle) and the E-CDFS 
Chandra area (large square). The radio and Spitzer areas used are 
all wider than the E-CDFS. The Herschel and the XMM-Newton 
areas are the boundaries of the "complete" and "broad" samples, 
respectively (see Sect. 13.31) . The sources of the "broad" sample 
are marked with blue symbols, whereas the sources with FIR 
flux density upper limits in the "complete" sample are marked 
with red symbols. 



luminosities and stellar masses (see Sects. [TT] and 13.2.1]) from 
the photometric redshift uncertainty is < 20%, therefore not 
important for the overall uncertaintie of the aforementioned 
values. 

The spatial limits of the difl'erent surveys described in this 
section are shown in Fig. [T] The grey- scale image is the com- 
bined 2-lOkeV image of XMM-Newton and the regions are 
the Herschel area limiting the combined GOODS-//^r^c/z^/-PEP 
catalogue (small rectangle), the 4 Ms Chandra area (the region 
where the efl'ective exposure is larger than half of its maximum 
value; solid circle), the XMM-Newton area (the region where the 
total integration time is higher than 1 Ms; dashed circle) and 
the E-CDFS Chandra area (large solid square). The radio and 
Spitzer areas described above are all wider than the E-CDFS. 
In this study, we use X-ray sources spanning the entire XMM- 
Newton region, and use information from all the other surveys 
to measure their star-formation rates and stellar masses. The 
smaller Herschel area is used to construct a "complete" sample 
of X-ray AGNs, where we have FIR detections or upper lim- 
its for the majority of the AGNs (see Sect. 13.3b , while for the 
wider area ("broad" sample) we use 70 yum measurements from 
the FIDEL survey. 



3. The sample 

In order to avoid high X-ray flux spurious detections due to rel- 
atively high background levels, we limit the sample to those 
sources which have a combined XMM-Newton exposure of 1 Ms 
or higher in the 2-10 keV band (356 sources in the main and 
supplementary catalogues of Ranalli et al. in preparation). To 
better constrain their positions we look for counterparts among 
the X-ray sources observed w ith the C handra surveys, namely in 
the 2 Ms CD FS catalogue of iLuo et al . (2010), the 4 Ms CDFS 
catalogue of 



iLehmer et al 



Xue et alJ ( Uoilh. and the E-CD FS catalogues of 



(12005') and ' ViranietaD (l2006h . The characteris- 
tic positional uncertainty of Chandra is < 1 arcsec, compared 
to the 4-5 arcsec of XMM-Newton. We keep Chandra counter- 
parts which are within 5 arcsec of the XMM-Newton position 
and find 311 unique associations. We also look for counter- 
parts in the 3.6 [im SIMPLE catalogue, using the likelihood ratio 
metho43 with a matching radius of 5 arcsec, and find another 
19 sources with LR > 0.85 and all with a rehability > 99.9^0. 
Most of them lie in the area not covered by the 4 Ms Chandra 
CDFS survey (see Fig. [T]). Our final X-ray catalogue contains 
330 sources with good positional constraints (< 1 arcsec) either 
from Chandra, or from Spitzer-\RAC. Of the 26 sources with no 
unambiguous Chandra, or Spitzer-\KAC counterpart, 23 are of 
low- significance (< 5 cr) and therefore likely spurious, and three 
are double sources in Chandra, not resolved by XMM-Newton, 
which we exclude from our sample. 

Next, we build the multi- wavelength catalogue of the X-ray 
sources using all the information available and the likelihood 
ratio method to select the counterparts, using the positional un- 
certainties provided in the various catalogues. We first combine 
the SIMPLE catalogue with both the ^-selected (Taylor e t^all 
2009) and the ^VT^-selected (iGawiser etal.Ll2006h MUSYC cat- 
alogues (preferring ^-selected sources in cases where they are 
detected in both catalogues) and find the optical-infrared coun- 
terparts of the X-ray sources, constraining their positions, and 
then we look for counterparts in the FIDEL and 24yum-prior 
Herschel catalogues. We find a counterpart in at least one of 
the optical or infrared catalogues for 328/330 X-ray sources; 
one source is too faint to be detected at any other wavelength 
than X-rays and the other is close to a bright optical-infrared 
source and is missed by the source detection algorithms. The 
positions of our opt ical-infrared coun terpa rts are in good ag ree- 
ment with those of lXue et al] (l201lh and iLuo etaL[(l2010h for 
the sources in common; more than 90% of the counterparts are 
within 0.7 arcsec of the optical positions given in those cata- 
logues. We i ise the optical po s itions to lo ok for radio counter - 
parts in the iKeUermann et al.l (l2008h and iMiUer et al.l (l2008h . 
and find 53 matches within a 2 arcsec radius, excluding X-ray 
sources which have multiple radio counterparts within 10 arcsec, 
which would be indicative of FRII radio AGNs. For the sub- 
mm catalogues, because of their large positional uncertainties 
(~ 8 arcsec) we use the likelihood ratio method to assign a 
FIDEL-24 counterpart to the sub-mm sources, and if this is the 



^ The likelihood ratio method ([Sutherland & SaunderslJl992h is usu- 
ally adopted in cases where a counterpart is sought in a crowded cata- 
logue (in this case the SIMPLE catalogue), and it uses the surface den- 
sity of objects of a given magnitude to estimate the probability that a 
counterpart at a certain distance is a chance encounter. An example of 
using this method to find infrared counterparts of XMM-Newton sources 
can be found in iRovilos et al.l (1201 1) 

^ The reliability is a measure of the probability that the selected coun- 
terpart is the correct one, and it is used in cases where more than one 
possible counterparts are found. 
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same as the EIDEL-24 counterpart of the X-ray source we con- 
sider it as reliable. We find a LAB OCA counterpart for five X- 
ray sources, and an AzTEC for two of these five. Finally, we look 
for redshifts and find 215 spectroscop ic re dshift determinations 
from the various catalogues listed in § 12.61 and 106 photometric 
redshift estimates. Nine sources have no redshift determination, 
because they are too faint at optical wavelengths. In the final cat- 
alogue we also include EIR upper limits for X-ray sources which 
are in the area observed by PACS (see Eig.[T]) with no detection. 
There are 155 XMM-Newton sources inside the RACS area and 
94 of them are detected. Of the remaining 61, 20 are in regions 
confused with nearby bright EIR sources, and 41 are upper lim- 
its. 

3.1. Stellar masses 

The most reliable method to derive stellar masses for galax- 
ies is the fitting of their broad-band spectral energy distri- 
butions (SEDs) with synthetic stellar templates with known 
star-formation histories and dust extinction properties (see 
IShaplev et all 1200 lUPapovich et al.L 12001. for the limitations of 
the method). The stellar component is important at optical wave- 
lengths (< 1 //m), but we use the full multi- wavelength informa- 
tion (excluding radio and X-rays) to fit the SEDs. The reason 
for this is that t he AGN can afi^ect the optical properties of the 
system (see e.g. iPierce et all [2010) , and by fitting a combina- 
tion of AGN and host templates using the infrared photometric 
information we can constrain the AGN contribution. 

Eor the opti cal-infra red SE D fitting we use the procedure 
described in iLusso et al.l (1201 ih . We appl y a minimisa- 
tion m ethod using stellar templates from the lBruzual & Charloi 
(l2003l) stellar synthesis code, applying solar metallicity, a mix- 
ture of constant and exponentially decaying star-formation rates, 
and a Galactic disk initial mass function (IME; Chabrier, 2003). 
We redden the stellar SEDs using the lCalzetti et al.l pOOO) law, 
and combin e the reddened SEDs with star-formation infrared 
SEDs from Chary & ElbazI (l200lh (105 templates with diff'er- 
ent E IR profiles in the 3 - 1000 yum range) and four AGN SEDs 
from ISilva et al.l (120041), which span from the optical to the 
far-infrared with diff'erent absorption properties (unabsorbed to 
A^H ^ 10^^ cm"^). Some characteristic results of the optical-to- 
infrared SED fitting, as well as a EIR-limit example can be seen 
in Eig. O These are examples of both AGN and starburst dom- 
inated SEDs. There is enough optical information (photometry 
and redshift) to fit an SED for 304 of the 330 sources. Sources 
with a detection in the lTaylor et al.l ('2009') catalogue (the major- 
ity of the optically-detected sources) typically have photometry 
in nine optica l - near-IR bandfl and those detected only in the 
iGawiser et al. I (12006) in six bands, which are used for the de- 
termination of the stellar masses. We do not take into account 
upper limits for the fitting, but check that the predicted flux den- 
sity value of the fitted SED is indeed lower than the limit. The 
reduced values of the best-fit models are typically in the 1-10 
range, after reprocessing the flux density errors using a quadratic 
combination with the 10%-level error, to account for the typical 
flux difl'erences between the SED templates used. 

The method we use to calculate the stellar masses induces 
uncertainties both from the choice of the difl'erent parameters 
fitted, and from the;^^ procedure itself. The derived stellar mass 
values are potentially strongly influenced by such uncertainties, 
especially at high redshif ts, like the maj ority of the sources in 
our sample (see also Michatowski et all 12012) . The uncertain- 
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Eig. 2: Examples of the SED fitting used to derive stellar masses 
and star-formation rates, selected to demonstrate the diversity in 
the SEDs of the X-ray sample. The star-formation component is 
plotted in red, the AGN component in blue, and the stellar com- 
ponent in green. In magenta we plot the combination fitted to the 
data-points. The stellar component is only confined by the opti- 
cal wavelengths where it usually dominates the flux, but there 
are cases where we detect a substantial AGN contribution. The 
star-formation component is confined by the EIR flux with rest- 
frame wavelength > 20 yum, and again there are cases where the 
EIR flux is dominated by the AGN. With downward arrows we 
plot far-IR upper limits (see Sect. 13.3b , where the resulting star- 
burst component is plotted with a dashed line. The numbering 
in the top-right corners of the panels refers to the preliminary 
XMM-Newton catalogue number. 



We do not use the photometry in the C/38 band for the SED fitting 



ties from the parameters included in t he stellar synthe s is pro - 
cedure are estimated to be ~ 0.15 dex (iBolzonella et al.L [2OIOI) . 
and checking the stellar masses and values of SED fits with 
diff'erent templates and diff'erent relative contributions, we esti- 
mate the final uncertainty in the stellar masses to be ~ 0.25 dex 
at the 90% confidence level (see iLusso et all 12012 '. for a more 
detailed description of the method). We also note that the use 
of the iChabrieil (12003 ) IME causes a s light und erestimation 
of the stellar masses with respect to the iKroupai (.2001) IME, 
in p articular they are on average lower by a factor of ~ 1.1 



(^seel Bolzonella et al.L[2QlQl: lPozzetti et al.'. '20 10: 'Mainline et all 
2011]). In this work we use the Chabrier (2003) IME in order 
to avoid an over-pre diction in the number of low-mass stars 
dHainline et al.L 1201 ll) . but we also combine the stellar masses 
with star-formation rates; the latter are based on infr^d lumi- 
nosities. This star-formation rat e proxy u ses the lKroupal (l200ll) 
IME for its calibration (see Murphy et al.Ll201 ll). so we increase 
the stellar masses we derive through the SED fitting by a factor 
of 1.1 to be consistent with the star-formation rates. 
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3.2. Star-formation rates 

Star formation in galaxies affects almost all of their observed 
properties, from the X-rays to the radio wavelengths, so there 
are traditionally a number of ways to measure the star-formation 
rate (SER). In the cases of AGN hosts we can rule out the X-rays, 
since they are completely outshone by the AGN (see also § 13.31) . 
In this work we test three methods based only on flux density 
measurements: i) infrared luminosity, ii) radio luminosity, and 
iii) optical SED fitting. 



3.2.1. Infrared luminosity 

The IR luminosity is arguably the most reliable tracer of star- 
forming activity and is well correlated with other tracers (see 
iKennicuttL Il998al: iKennicutt & EvansL I2Q12L for reviews). The 
IR photons are emitted by the dust surrounding young stars, 
which is heated by their ultra-violet radiation. In this paper we 
will use the integrated rest- frame 8 - lOOOyum luminosity and 
the equation: 



I SER 
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;x 10 



-44 



erg s" 
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from lMurphy et al] (l2Qlll) . In order to measure the IR luminosity 
we perform the SED decomposition described in § 13.1! anew, us- 
ing only t he infrared data - points from Spitzer and Herschel, the 
complete IChary & ElbazI (1200 ih host templates (i.e. including 
wavelengths < 3 //m), the AGN templates, and ignoring the syn- 
th etic stellar part. The me t hod we use is the same as described 
in iGeorgantopo ulos et al.l ('2011a') and 'Georgantopoulo s et al.l 
([2011b)," using the SED templates described in § 13.11 We do this 
in order to avoid the degeneracies in the optical wavelengths 
between the stellar and AGN light, which could affect the fit- 
ted AGN contribution to the infrared luminosity; this way we 
fit fewer free parameters. Eor the infrared SED decompositions 
we require at least three mid-IR points from Spitzer-IKAC to 
determine the shape of the mid-IR part of the SED, and at least 
one flux density determination in a rest- frame wavelength higher 
than 20 yum, which we can use to constrain the far-IR part of the 
SED; there are 125 X-ray sources which comply with these cri- 
teria. The EIR flux density determination comes from EIDEL- 
70yum (for z < 2.5), R\CS-100yum (for z < 4), RACS-160yum, 
sub-mm (LABOCA and/or AzTEC), or a combination of them. 
We use the galaxy component (AGN-free) to calculate the star- 
formation rates, and for 14 cases the AGN component dominates 
even the longest wavelength IR data-point available, so the de- 
termination of the AGN-free part of the IR emission is not reli- 
able. The uncertainties in the IR luminosity values come mostly 
from the SED decomposition, and a check of the values of 
fitting secondary solutions not selected yields an uncertainty of 
-0.3 dex (or a factor of two) in the 90% confidence level, for 
the vast majority of the sources. The uncertainties arising from 
the far-IR flux errors are much lower. 



3.2.2. Optical SED fitting 

The star-formation rate can be derived as a by-product of 
the optical SED fitting performed in § 13.11 using the star- 
formation history and age assumed, and the normalisation from 
the photometry. Similar methods ha ve been widely used in deep 
fields, including the CDES (e.g. Brusa et aU i2009i) . especially 
if far-infrared photometry is not available. In the next para- 
graph we will test its reliability, since it is a highly model- 
dependent method with systematic uncertainties arising mainly 



fro m the IME, star-formati on history and extinction law used 
(^see'Bolzonella et al.L [2QTQb . 

In Eig.[3alwe plot the star-formation rate measured from the 
infrared luminosity of the sources detected in the far-infrared 
against the SER measured from the optical SED fitting after cor- 
recting for extinction, excluding the AGN contribution for both 
cases. All X-ray sources with an infrared measurement with rest- 
frame wavelength above 20 jim and an optical identification (in 
MUSYC) are plotted with a circle. We exclude these 14 cases 
where, according to the SED decomposition, the flux density 
of the longest wavelength data-point is dominated by the AGN, 
so that the SER cannot be constrained (see Eig. |2c]). In Eig. [3a| 
we also include SER upper limits for the X-ray sources in the 
Herschel area which are not detected by Herschel\ for one of 
the 41 X-ray sources with EIR upper limits we do not have any 
photometric data-points to perform an SED fitting in the mid- 
infrared, while another four are not detected in the optical, or 
have no redshift determination, so 36 upper limits are plotted in 
Eig.[3al The solid line is the 1:1 line and the dotted lines mark 
the ±1 dex region. We can see that out of the 109 points of Eig. 
[3al 79 are between the dotted lines, while for 29 sources using 
the optical SEDs underestimates the SER by more than an order 
of magnitude; for one (and two upper limits) the SER is over- 
estimated by more than an order of magnitude, assuming that 
the infrared SER is reliable. The star-formation rate estimated 
from the optical SED is a highly model-dependent value, and is 
very sensitive to the star-formation history assumed in the stellar 
synthesis models, which shape the optical SED. It is also sen- 
sitive to dust extinction, which would cause an underestimation 
of the SER, explaining the behaviour we see in Eig. [3al Due to 
this large scatter and systematic offset, we do not rely on the op- 
tical SED fitting to derive star-formation rates and use it only 
for stellar mass determinations. The stellar mass is an integrated 
value and therefore less sensitive to the assumed star-formation 
history. Moreover, we do not detect any obvious dependence of 
the difference between the SER determination using the two es- 
timators on X-ray or optical classes, which would indicate AGN 
contamination as the cause of the scatter. 



3.2.3. Radio luminosity 

The radio luminosity of star-forming g alaxies is tight ly corre- 
lated with their infrared luminosity (see Condon, 1992, for a re- 
view), and this correl ation holds even for cosmologically signif - 
icant redshifts (z ~ 2: lAppleton et al.Ll2004l:llvison et al.[l2010h : 
we test it in this work as a possibility to derive the star-formation 
rates in X-ray sources without far-infrared detections. The radio 
emission in star-forming systems is generated by synchrotron 
radiation from relativistic electrons accelerated by supernova in- 
duced shocks and free-free emission in H ii regions. The caveat 
is that AGNs themselves can produce radio emission through 
radio jets or compact high surface-brightness synchrotron core 
emission from relativistic electrons heated by the AGN. There i s 
a dichotomy in the radio power of quasars (Miller et al.L Il990h . 
with sources having Lsghz ^ 10^^ WHz"^ being characterised 
as "radio-loud" and their power source being closely connected 
to the AGN, and sources having Lsghz ^ lO^^WHz"^ being 
characterised as "radio-quiet" and having a controversy about 
their power source. Alternatively, the radio to optical flux ra- 
tio is used in some stu dies to differentiat e between radio-quiet 
and radio-loud AGNs (iKellermann et al.L Il989 ). However, the 
dichotomy between radio-loud and radio- quiet sources is no t 
clear if more complete samples are used (^e.g. lWhite et al.Ll2000h . 
with many objects in the "intermediate" region, making the tran- 
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Fig. 3: Comparison of the star formation rates of X-ray sources, measured using the infrared 8 - 1000 jum luminosity, and the optical 
SED fitting (a), or the 1.4 GHz radio luminosity (b). Black circles mark reliable fits in both optical and far-infrared (or radio) 
wavelengths, excluding cases where the AGN SED dominates even the longest wavelength flux density available. Sources spanning 
over the whole 3 Ms XMM-Newton region are used in these plots, while the downward arrows refer to Herschel non-detections in 
the PACS area (35 cases), and leftward arrows in the right panel refer to radio upper limits spanning over the whole XMM-Newton 
region (112 cases). The solid and dashed lines in both panels represent the 1:1 relation and the ±1 dex deviation, respectively, and 
the vertical dashed line in the right panel marks the calculated star-formation rate of a source having Li.4ghz = 10^^-^ WHz"^ = 
2Q31.5 gj,g g-i The significant deviation form the 1:1 relation of both cases and the significant scatter of the optical SED (left) 
case prevent us from using any of the two SFR tracers in this paper. 



sition smooth with only a vague limit. Recently. iPadovani et al.l 
(l2QTTh . using the luminosity functions of diff'erent types of ra- 
dio sources in the CDFS, argue that the major contribution 
to radio power in the diff'use part of radio-quiet AGNs comes 
from star-formation, though taking a somewhat stringent limit to 
characterise the radio sources based on their radio luminosities 
(^i.4GHz = 1024-5 WHz-i = 10^^-5 erg s-iHz-\ combined with 
other observational characteristics). 

For this work we test how reliable the radio luminosities 
are in estimating star-formation rates for an AGN sampl e using 
the VLA 1.4 GH z flux densities from iKellermann et al.l (l2008l) 
and Miller et al. ( 2008), also checking the VLBI catalogue of 
iMiddelberg et al. ( 2011) to exclude any high surface-brightness 
compact cores, characteristic of non-thermal nuclear emission, 
not c onnected to the star formation (e.g. Giroletti & Panessa 
[2009h . We calculate the radio luminosities using 



according to the Li.4ghz > 10^^-^ erg s ^ Hz 



^1.4GHz = 47r<ii S l.4GHz(l + z)' 



a-l 



(2) 



where a is the radio spectral index, assuming Sy oc y"^, and 
it is calculated from the relative radio flux densities at 1 .4 and 
5 GHz. In cases where the 5 GHz flux density is not available 
we assume a = 0.8, characteristic of synchrotron emission (see 
ICondonLfT992h . There are 53 X-ray sources with a radio counter- 
part within 2 arcsec and without another radio source closer than 
lOarcsec, the latter would suggest an FRII radio-loud source. 
Eight of these sources have a high surface-brightness core de- 
tected with VLBI with a flux density above 0.5 mJy, and are re- 
moved from the test sample, and a further eight are radio-loud 



criterion (three are 

also detected in the far-infrared and are included in Fig. [30). The 
star-formation rate is calculated using 



SFR 

Mo yr-i 



6.35 X 10"^^ 



-1.4 GHz 



erg s" 



Hz- 



(3) 



(IMurphv et al.L l201lh . In Fig. [3b] we plot the star-formation 
rates from the infrared and radio luminosities for sources be- 
ing detected in both bands, keeping the same range and sym- 
bols as in Fig. [3al We also plot the radio and infrared (for 
the Herschel area) upper limits with arrows. The scatter in this 
case is significantly lower than in Fig. [3a| However, the mean 
log(SFRiR/SFRi.4GHz) of the radio detections is -0.48 (with- 
out taking into account the VLBI sources and the upper lim- 
its) with a standard deviation of cr = 0.39. The dashed verti- 
cal line in Fig. [3b] marks the calculated SFR of a source hav- 
ing Li.4GHz = 10^^-^ erg s"^ Hz~\ thus being border-line radio- 
loud according to the limit of 'Padovani e t al.l (1201 1 ). If we keep 
this limit and calculate the mean log(SFRiR/SFRi.4GHz) of only 
the radio-quite objects, its mean and standard deviation become 
-0.45 ± 0.40, so there is still contamination from the AGN emis- 
sion; most probably we are detecting in the radio band those 
sources which are in the top of the radio flux distribution, some- 
thing which is supported by the location of the upper limits. 
Because of this contamination, we do not use the radio power 
as a star-formation proxy in our AGN sample. 
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Fig. 4: Optical magnitude against 2-lOkeV flux density for the 
XMM-Newton sources with robust SFR estimations. The lines 
mark the -1 < log(/x//opt) < 1 region and open symbols sources 
with Lx < 10"^^ erg s"^ We exclude the 10 sources which have 
both Lx < 10^^ erg s"^ and log(/x//opt) < -1 from our final sam- 
ple, as normal galaxy candidates. 
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Fig. 5: Histograms of the basic properties (star-formation rates, 
stellar masses, specific star-formation rates, and redshifts) for 
the 99 hosts of the X-ray AGN sample used in this work. The 
shaded histograms show the same properties for X-ray QSOs, 
i.e. sources with Lx > 10"^"^ erg s"^. 



3.3. Final sample 

We start with a sample of 356 X-ray 2-lOkeV selected sources 
from the 3 Ms XMM-Newton survey with a total integration time 
more than 1 Ms, 330 of which have good (< 1 arcsec) positional 
constraints from Chandra or Spitzer. We have enough optical 
information to fit an SED and calculate the stellar masses for 
304 of these 330 sources. On the other hand, we rely on the 
infrared flux density to constrain the SFR of our sample, and 
there are 111 sources with an SFR measurement from the FIR 
flux. For 109 of them we can also calculate the stellar mass from 
the optical SED. 

As we are dealing with faint X-ray fluxes, there is the possi- 
bility that for some of the X-ray sources the X-rays trace normal 
star-forming galax ies instead of the AGNs (i Ranalli et al.L 120031: 
iBauer et al.U2004l) . In Fig.|4]we plot the hard X-ray flux densities 
against their optical (7?-band) magnitudes. The lines mark the 
-1 < log(/x//o pt) < 1 region where the bulk of the AGNs are ex- 
pecte d (see e.g. lStocke e'tailll99ll:lElvis et al.L [19941: IXue et all 
l201lh . Sources with log(/x//opt) < -1 are candidates for being 
normal galaxies instead of AGNs (see iTzanavaris et al.l 120061: 
^Georgakakis et al.L l2006l) . Moreover, most normal galaxies have 
X-ray luminosities not exceeding 10^^ erg s"\ except for a few 
extremely star-formi ng sources, mainly detected in sub-mm 
wavelengths (see e.g. [Alexander et al.L I2QQ5I: iLaird et ani2010h . 
Sources with luminosities below the 10"^^ erg s~^ limit are plotted 
with open circles in Fig. (4] There are 10 sources compliant with 
both the /x//opt and the Lx criteria, and they do not show any 
signs of obscuration in their X-ray spectra, so we remove them 
from the AGN sample. 

A fundamental property of each galaxy is its specific star- 
formation rate (sSFR), which is defined as the ratio of its star- 
formation rate to its stellar mass. It is indicative of how efficient 



the galaxy is forming stars. To calculate the sSFR for the AGNs 
in our sample, we use the star-formation rates measured from the 
infrared luminosity. We have calculated the sSFR for 99 X-ray 
AGNs, 77 with spectroscopic redshift and 22 with photometric 
redshift. The SFR, stellar mass, sSFR and redshift histograms 
are shown in Fig. [5] 

The sample described in the previous paragraph (hereafter 
the "broad" sample) is not a complete sample of X-ray detected 
AGNs, because of the various selections of the sources, which 
limit their number from 356 to 99. This incompleteness might 
aff'ect the statistical properties. In order to account for that, we 
create a more complete sub-sample constrained in the Herschel 
area where we have the most sensitive far-IR measurements. In 
this area (marked with the small rectangle in Fig.[T]) for the X-ray 
sources for which a far-IR counterpart is not found, a FIR up- 
per limit is calculated from the sensitivity map of the GOODS- 
Herschel-PEP survey. Thus, we have 155 X-ray sources, 94 of 
them are detected in the FIR, and for 41 we can calculate an 
upper limit to their lOOyum and 160yum fluxes. Twenty sources 
lie within a 10 arcsec region of a nearby bright FIR source and 
an upper limit cannot be calculated, they are however a ran- 
dom sub-sample, not afl'ecting the completeness. Out of the 135 
(155-20) sources, eight are associated with normal galaxies not 
hosting an AGN, a further 12 do not have sufficient optical or 
mid-infrared information, or a redshift estimate for an SED fit, 
and for a further seven the emission from the AGN dominates 
over the FIR flux (the AGN-related flux go the highest wave- 
length data-point is higher than the star-formation related), mak- 
ing a star-formation measurement not reliable. Summing up, we 
calculated the SFRs and stellar masses of 108 out of the 127 
(^ 85%) X-ray AGNs in the GOODS-Herschel region, for which 
a FIR flux determination is possible (76 detections and 32 up- 
per limits). Hereafter we will call this the "complete" sample. 
The basic properties (2-10 keV luminosities and redshifts) of the 
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Fig. 6: 2-10 keV luminosities against redshift for all the X- 
ray sources with redshift determinations in the 3 Ms XMM- 
Newton survey. The original parent sample is plotted in grey 
crosses, while black and red filled circles are plotted for the 
X-ray AGNs with both a stellar mass and a SFR measure- 
ment ("broad" sample). Red symbols represent sources in the 
GOODS-//^r^c/z^/-PEP area, filled for Herschel detections and 
open for FIR upper limits ("complete" sample). The line is the 
6.5 X 10"^^ erg s"^ cm"^ flux limit, assuming F = 1.7. 

complete sample are shown in Fig. [6] with red symbols (filled and 
open circles for FIR detections and limits, respectively), while 
the properties of the overall sample are plotted in black symbols, 
and the rest of the X-ray sources are plotted in grey crosses. 

4. Results 

4.1. sSFR-U 

In previous studies there has been a controversy about the exis- 
tence of an observational connection between the AGN and the 
host galaxy activity. In Fig. [7] we plot the SFR against the hard 
X-ray luminosity of the 99 X-ray AGNs with an estimate of the 
SFR described in the previous section and the 32 upper limits. 
There are also two X-ray sources with lower limits in their X- 
ray luminosities. These are Compton-thick sources whose X-ray 
spectra are dominated by a reflection component according to 
the spectral fits, and their unobscured luminosities cannot be de- 
termined. The lower limits in Fig [8] are their observed luminosi- 
ties. Because of the limits, for the statistical anal ysis we use the 
ASURV package (Rev. 1.3 LaVall ev et al.Lll992h. which im ple- 
men ts the methods presented in Feigelson & NelsonI (11985 ) and 
llsob e et al. (1986). Using the generalised Kendall's r method in 
order to include upper limits, and all the data-points of Fig. [71 
we find that the SFR is strongly correlated with the hard X-ray 
luminosity, with a null hypothesis probability lower than 0.01%. 
To simulate a mass-matched sample and study the activity of the 
host independent of its size, we calculate the specific SFRs of 
the sample and plot it against the X-ray luminosity in Figure O 
Performing the same method, we find again that the two values 
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Fig. 7: Star-formation rate against hard X-ray luminosity for 
X-ray selected AGNs. We plot values of 99 far-infrared de- 
tected X-ray AGNs, and 32 FIR upper limits in the area cov- 
ered by the deep Herschel survey. The dotted line is the expected 
FIR luminosity of a p ure- AGN source, translated into SFR (see 
'Mullanev et alj,|2011|). We find a strong correlation between the 
star-formation rate and the X-ray luminosity. 

are strongly correlated. However, iMullaney et al.l (l2012ah have 
shown that the correlation between the X-ray and the infrared lu- 
minosity is sensitive on the evolution of the infrared luminosity 
with redshift, and this might be aff'ecting the sSFR-Lx correla- 
tion we observe here. To test this hypothesis we apply the partial 
correlation test of Akritas & Siebert (1996) and find significant 
correlations of the sSFR with both redshift and X-ray luminosi- 
ties at the 4.6cr and A.Acr levels, respectively. However, this par- 
tial correlation test tends to incorrectly reject the null-hypothesis 
in cases where the two "independent" param eters (here a nd 
z) are also correlated with each other (see Ke lly et aP. l2007l) as 
in this case, which limits the reliability of the test. 

In order to further check the eff'ect of the redshift in the 
sSFR-Lx correlation, we divide our sample into six redshift bins, 
containing a roughly equal number of data-points (21 or 22), and 
repeat the Kendall's r method for each bin separately and for 
different combinations. The description of the redshift bins and 
the results of the test (null hypothesis probability) are shown in 
Table[TJ We can see that there is no correlation at lower redshifts, 
but there is a possible correlation for z > 1.15 (bins 4, 5, 6) with 
> 95% significance. If we merge adjacent bins in order to in- 
crease the number of data-points in each bin, the correlation is 
again found for z > 1.12 (bins 4-5, 5-6) with > 99.8% signifi- 
cance. In this case however the redshift dependence is not negli- 
gible. We also note that within the redshift bins the sample is al- 
most luminosity-limited. To further check if the X-ray flux limit 
afl'ects the previous result, we exclude sources with X-ray lumi- 
nosities lower then the luminosity limit of the highest redshift 
limit of each redshift bin, to create truly luminosity-limited sub- 
samples. Repeating the analysis in those sub-samples, the pre- 
vious result does not change significantly, except in the highest- 
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Fig. 8: Specific star-formation rate against hard X-ray luminosity 
for X-ray selected AGNs. The grey, blue, green and red symbols 
refer to z < 1.120, 1.120 < z < 1.615, 1.615 < z < 2.455 and 
z > 2.455, respectively. The error-bars and their associated lines 
refer to the mean luminosities and specific SFRs of the high- and 
low-luminosity bins within each redshift bin, using the Kaplan- 
Meier estimator (see Sect. 14.1b . We do not detect a significant 
correlation between the X-ray luminosity and sSFR for the low- 
est redshift bin, but do detect a significant correlation for higher 
redshifts (z> 1). 



redshift bin (z > 2.305). In Fig. [8] we colour-code the data-points 
with respect to their redshifts, in grey we plot bins 1-2-3, in blue 
bin 4, in green bin 5 and in red bins 6; the results of this binning 
are shown in the last column of Table [T] We divide each data 
compilation into a low-luminosity and a high-luminosity bin in- 
cluding an equal number of sources, and plot the mean sSFR 
and its associated error, calculated using the Kaplan-Meier esti- 
mator in ASURV, and the mean luminosity of the bin. The sSFR- 
Lx correlation is evident for the blue, green and red data-points 
(z> 1.15). 

A possibly important factor aff'ecting the previous analysis 
is the FIR selection of the sources of our final "broad" sam- 
ple, which reduces the number of X-ray sources from 356 to 
99, being biased in favour of sources with higher SFRs. In or- 
der to check whether this has an eff'ect on the apparent sSFR- 
Lx correlation, we repeat the previous analysis in the small area 
covered by the GOODS-Herschel survey (the "complete sam- 
ple"; see Figs.[T]and[6l). In this case, the analysis is performed in 
broader redshift bins due to the smaller number of sources, and 
the results are similar to those described in the previous para- 
graph; there is no sign of a sSFR-Lx correlation below z < 1.2, 
but above this redshift the correlation is > 95% significant. 

We note here that there is a small number of X-ray sources, 
which are detected in the far-infrared in a rest- frame wavelength 
A > 20 fim, but its flux density is dominated by the AGN emis- 
sion, according to the SED decomposition performed. There are 
14 such cases in the "broad" sample and seven in the "complete" 
sample with redshifts z ~ 1-3. These sources could populate 
the low-(s)SFR - high-Lx area, however their far-IR luminosi- 



ties cannot be constrained, not even with an upper limit. If we 
consider the SFR estimations from the optical SEDs, although 
unreliable, they are consistent with the (s)SFR-Lx correlation. 
Moreover, their number is ~ 10% of the sample used, so we are 
confident that they will not aff'ect the result. 

4.2. Redshift evolution 

The average sSFR of star-form ing galaxies increases w i th red - 
shift at least up to z ~ 2 (Elbaz et al.L l2007l: iDaddi et al.1 l2007h . 
and in this section we investigate how the hosts of an AGN 
evolve with respect to the general population. In Fig.[9a|we plot 
the sSFR of the AGN hosts against the redshift. The vertical lines 
refer to the redshift bins of Table [T] while the solid curve is the 
expected main-sequence sSFR, according to 



sSFRMs[Gyr-i] 



26xr2-2. , 

cosmic ' 

2, 



z< 2.156 
otherwise 



(4) 



where ^cosmic is given in Gyr. The grey area denotes the borders of 
the starburst and quiescent areas, de fined as double an d half the 
main sequence sSFR, respectively (lElbaz et al.Ll201ll) . We note 
here that the increase of the main- sequence sSFR does not con- 
tinue forever, and in Elbaz et al. ( 201 1) the density of data-points 
supporting the above relation dramatically decreases at z > 2.5. 
There is evidence that the main- sequence sSFR i s constant above 
z ~ 2 (IStark et al.L 1200 91: 'Gonzalez et al.", '2010) with a value of 
sSFRms ~ 2 Gyr~^ According to the above relation the value of 
sSFRms = 2 is reached at z = 2.156, so above this redshift we 
assume a constant relation. There is a hint that the hosts of the 
AGNs in our sample are mostly in the main-sequence and star- 
burst regions, while they generally follow the behaviour of the 
main sequence with redshift. With green crosses we mark the po- 
sitions of X-ray QSOs having intrinsic L2-iokeV > 10"^"^ erg s"^. 
The majority of them (19/25 sources) are consistent with being 
in the starburst region with sSFR/ sSFRms (^) > 2. 

This behaviour is also evident if we plot the "starburstinesfl" 
against redshift in Fig.[9bl The "starburstiness" is the ratio of the 
sSFR of the source over the main-sequence value at the given 
redshift. The vertical dashed lines in Fig.|9b]are identical to those 
in Fig. [9al and the grey area again marks the main sequence. 
The symbols of the data-points are identical to Fig. |9a| (black 
for all the sources and green for X-ray QSOs). In each redshift 
bin we also show the average "starburstiness" and its associated 
statistical uncertainty calculated using the Kaplan-Meier estima- 
tor. For the QSO case we have re-binned the data into four red- 
shift bins to improve the statistics of the sample. The behaviour 
seems to depend on redshift: the QSOs in the first redshift bin 
with 0.976 < z < 1.499 have an average sSFR consistent with 
that of the overall population, while higher redshift QSOs are on 
average more "starbursty", having sSFR more than double that 
of the main sequence. 

To check how much the complex source selection aff'ects 
those results, we repeat the previous analysis for the "complete" 
sample, where we have FIR upper limits for most of the X-ray 
sources. We use the same technique, and the result is shown 
with the dashed lines in Fig. |9bl It is consistent within the sta- 
tistical uncertainty with that of the "broad" sample. For the QSO 
case however, the diff'erence between the complete and the broad 
sample is significant in the first redshift bin, owing to the scarcity 
of such objects. The mean sSFR of the QSOs seems to be in the 
main sequence for z < 2 and in the starburst region for higher 
redshifts. 



starburstiness = sSFR/ sSFRms (z) 
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Table 1: Results of Kendall's r method for the correlation between the specific star-formation rate and the hard X-ray luminosity 
for diff'erent redshift bins. The null hypothesis probability in each redshift bin is shown in column 3 and in combinations of bins in 
columns 4 and 5, with a lower null hypothesis probability meaning a tighter correlation. 



Bin number of sources redshift range Null Hypothesis (%) Null Hypothesis (%) Null Hypothesis (%) 



21 
22 
22 
22 
22 
22 



0.000 - 
0.625 - 
0.759 - 
1.156- 
1.605 - 



0.620 
0.755 
1.113 
1.599 
2.299 



> 2.305 



34 

8.5 

43 

5.4 

0.97 

5.2 



7.9 (bins 1-2) 
4.3 (bins 2-3) 
9.5 (bins 3-4) 
0.18 (bins 4-5) 
0.11 (bins 5-6) 



1.5 (bins 1-3) 

5.4 (bin 4) 
0.97 (bin5) 
5.2 (bin 6) 
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Fig. 9: Evolution of the specific SER and the "starburstiness" with redshift for X-ray selected AGNs with long wavelength informa- 
tion. The dashed vertical lines refer to the redshift bins of Table [T] the solid line in panel (a) is the expected main sequence sSER 
according to Equation |4l and the grey area s denote the limits of the starburst and quiescent areas. Their borders are double and 
half the main-sequence sSER, according to Elbaz et al Green crosses mark the positions of X-ray QSOs having intrinsic 

^2-iokeV > 10"^"^ erg s"^. In panel (b), the solid lines and respective data-points are the running means of the "starburstiness" for X- 
ray AGNs and QSOs in the "broad" AGN sample, and the dashed lines are the running means of the "complete" X-ray AGN sample 
(see Sect. 13.31) . There is a general trend for the sSER of AGN hosts to follow the main sequence, so that the median "starburstiness" 
is constant with redshift. The QSO hosts on the other hand have sSERs which are somewhat higher. 



In order to further investigate how the "redshift eff'ect" af- 
fects the correlation found between the sSER and the hard X-ray 
luminosity, in Sect. 14. II we plot the starburstiness defined in this 
section against the X-ray luminosity in Eig. \T0\ The results of 
the statistical analysis (null hypothesis probability) of the same 
redshift bins as in Table [T] are presented in Table O The analy- 
sis shows again a > 95% correlation for redshifts z > I and no 
correlation at lower redshifts. We therefore assume that the cor- 
relation between the host and galaxy activity is not aff'ected by 
the evolution of the infrared luminosity with redshift. 



4.3. sSFR-Nu 

In earlier studies there have been some hints of a correlation 
between the star-forming activity of the h ost and the AGN ob- 
scuration in the X-rays. iPage et al.l (l2004l) found that X-ray ab- 



sorbed sources are more likely to be detected at sub-mm wave- 
lengths because of their extreme star-formation rates, although 
the abs orbed AGN sample consisted only of type I (broad-line) 
QSOs (Pageetal., 2001), which are not a representative sam- 
ple. Moreover, Alexander et al. ( 2005) found that the majority 
of radio-detected SCUBA sub-mm sources are consistent with 
being heavily obscured AGNs, with Nu ^ 10^^ cm ~^, although 
the active nucleus is not bolometrically dominant. iBauer et aP 
( 2002) found hints that X-ray sources with sub-mJy radio coun- 
terparts (tracing star formation) are, on avera ge, more obscured 
than the unmatched po pulation, confirmed by iGeorgakakis et al.l 
(i200 4). Subsequently, Rovilos et al.l (l2007 l) using a combination 
of the 1 Ms CDES and the E-CDES surveys found that such a 
trend was confined only to AGNs with any evidence for X-ray 
obscuration (A^h > 10^^ cm"^), linking it with line-of-sight ef- 
fects. However, deeper surveys both in X-rays and at infrared 
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Table 2: Results of Kendall's r method for the correlation between the starburstiness and the hard X-ray luminosity for different 
redshift bins. The null hypothesis probability in each redshift bin is shown in column 3 and in combinations of bins in columns 4 
and 5. 



Bin number of sources redshift range Null Hypothesis (%) Null Hypothesis (%) Null Hypothesis (%) 



21 
22 
22 
22 
22 
22 



0.000 - 
0.625 - 
0.759 - 
1.156- 
1.605 - 



0.620 
0.755 
1.113 
1.599 
2.999 



> 2.305 



84 
6.7 
69 
6.0 
1.0 
5.2 



31 (bins 1-2) 
4.5 (bins 2-3) 
29 (bins 3-4) 
0.5 (bins 4-5) 
3.7 (bins 5-6) 



4.4 (bins 1-3) 

6.0 (bin 4) 
1.0 (bin 5) 
5.2 (bin 6) 



o 



- 1 1 1 1 M 1 1 
• 

r I 


1 1 1 1 1 M 1 1 1 1 1 1 1 
• • 


M 1 1 1 1 1 1 1 1 L 
• 

•• • 






— • 


• |_ 






• 

.4, 

1 1 1 1 1 1 1 1 


nI 

^ • • 

1 1 1 1 1 1 1 1 1 


• 

1 1 1 


_i_i 


K • : 

III 1 1 1 1 1 1 1 1 



d 10*2 10^3 10^4 10^5 



logLg-iokev (ergs-1) 

Fig. 10: "Starburstiness" against (2-10) keV X-ray luminosity 
for the sources in our sample. The grey area is the same as in 
Fig. [9] and the data-points and error-bars are as in Fig. [8l sub- 
stituting the gray point with cyan. We again do not detect any 
significant correlation between the "starburstiness" and the X- 
ray luminosity for z < 1, but there is a correlation for z > 1, 
within the redshift bins. 



wave l engths failed t o reproduce those results ( e.g. iLutz et all 
120 Id: iRosario et al.l l2012l: iTrichas et al.L l2012h . Here, we use 
the deepest XMM-Newton survey, providing good quality X-ray 
spectra, combined with the deepest Herschel R\CS observa- 
tions, and a n SE P decomposition technique to clarify this is- 
sue. In Fig. Illal we plot the hydrogen column density of the 
sources for which we have good X-ray spectra against their star- 
formation rates. We apply a value of A^h = 10^^ cm"^ to sources 
which show no signs of obscuration and A^h = 5 x 10^"^ cm"^ 
to Comp ton-thick AGNs. We do not find any significant correla- 
tion between the two values. Using the Kendall's r method we 
find a null-hypothesis probability of 90%. Excluding unobscured 
AGNs or splitting the data into redshift bins does not change 
this result; the null hypothesis probability is always higher than 
20%. Neither is there a significant correlation in the high lumi- 
nosity AGN (Lx > lO'^'^ergs"^) or the high redshift (z > 1.5) 
sub-samples. In order to simulate a mass-matched sample and 
correct for any redshift eff ^ects in the column density and sSFR 
values (e.g. iHasingeil 12008.) we also plot the specific SFR and 



the starburstiness against A^h in Figs. Illbl and lllci Performing 
all the previous tests, we again do not find any significant cor- 
relation, except for the sSFR and starburstiness of AGNs with 
0.7 < z < 1.4, where we find hints of an anti-correlation 
at the 95% level. However, considering the behaviour of the 
overall sample and the complex selection eff'ects to shape that 
sub-sample, we do not consider it important. This behaviour is 
in broad agreement with models assuming a cl umpy absorber 
(Elitzur & ShlosmanL [ 20061: iNenkova et al.Ll200"8b . where the ab- 
sorption strongly depends on the number of absorbing clumps 
crossing the line-of-sight, however there is a number of obscured 
AGNs (with A^H > lO^^-^ cm'^) and sSFR/sSFRMsfe) < 1 which 
are still hard to explain with these models. 

4.4. Rest-frame colours 

The colour-magnitude diagram (CMD; the rest- frame U - 
B or U - V colour plotted against the absolute B or 
V magnitude) is used in a number of studies to check 
the evolutionary stage of the AGN hosts. The hosts of 
AGN s are concentrated in and around the "green val- 
ley" (Nandra et al., 2007; Rovilos & GeorgantoDoulos, 20071; 
Silver man et al., 2008; Georgakakis et al., 2008; Hickox et al.,, 
20091: iGeorgakakis & Nandral l2Qlll) . which is thought to sign- 
post the transition phase from a starburst to a "dead" ellipti- 
cal. There are however a number of factors that aff'ect the po- 
sition of a source (especially an AGN) in the CMD making 
the meaning of the above observation unclear. For example, it 
has been noted that both the AGN contribution and dust ob- 
scurat ion can alter the observed optical colours of the AGN 
hosts (iPierce et al.Ll201Ql:ICmiamone et alll2QlQbl:lLusso et al.L 
HoO), making them bluer or redder. Moreover, AGNs are usu- 
ally found in relatively high stell ar mass hosts (M^ ^ 10^^ - 
IO^Mqi Kauff'mann et al, 2003|; iBrusa et al.L I2QQ9I: IXue et al.L 
l201Ql:lMullaney et al.i. .2Q12a. this study), and that in turn makes 
them "avoid" the blue cloud; it is observed that the concentra- 
tion of AGNs around the "g reen valley" is not detected when 
using mass-matched s amples (ISilverman et al.l I2QQ9I: IXue et al.l 
l2QlQt iMuUanev et al.L [2QT2ah . In this section we use the sam- 
ple of AGNs for which we have an independent way to measure 
the star-formation activity, to check the validity of the colour- 
magnitude diagram without correcting the optical magnitudes 
for the host contribution or dust reddening. In addition, we use a 
colour-mass diagram instead of the colour-magnitude approach, 
in order to simulate a mass-matched sample. 

The positions of the red sequence, blue cloud and green 
valley in the colour-magnitude (and colour-mass) diagram are 
strongly dependent on redshift (iBellL l2QQ3l: iBorch et al.L l2006t 
Peng et al., 2010). In this section we use the U - B rest- frame 
colours of a sample of AGN hosts spanning from z ~ 0.1 to 
z ~ 4 and we use the parametrisation of Peng et al.l (l201Qh to de- 
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Fig. 11: Star-formation rate, specific star-formation rate, and 
"starburstiness" plotted against the hydrogen column density for 
65 sources with XMM-Newton spectra and robust SED fitting 
(including nine with Herschel upper limits). The grey area in 
panel (c) is the main-sequence, as in Figs.|9]and[T0l The errors 
on the hydrogen column density are in the range of < 10% to 
30%. We do not detect any significant correlation between the 
two values, even if we split the sample into redshift or luminos- 
ity bins, in any of the three plots. 



fine the dividing li ne between blue an d red galaxies, extrapolated 
to higher redshifts.lXue et al.) (l20ld ) have demonstrated that the 
colour bi-modality of non-AGNs exists up to at least z ~ 3. We 
use the fitted SEDs including the AGN contribution and the fil- 
ter curves of the COMBO- 17 survey (Bell, 2003) to measure the 
optical colours of the AGN hosts. In Fig. [12] we plot the star- 
burstiness of the AGN hosts against their "redness", defined as 
the deviation of their rest-frame colours from the dividing line 
in the colour-mass diagram. We observe a clear anti-correlation 
between the two values, which is statistically significant at a 
> 99.9% level, ^ ye parametrise it using t he Buckley- James re- 
gression method (iBuckley & Jamesll 19791) : 




log 



sSFR 



-0.5 0.5 

**redness** = (U-B)-(U-B)Rs 

Fig. 12: Starburstiness against deviation from the dividing line 
between red and blue galaxies for the sources of our sample. 
With blue and red symbols are plotted galaxies in the blue cloud 
and the green valley, and red- sequence galaxies respectively; 
the dividing line (the border of the red sequence) is the d ashed 
vertical line, calculated from the parametrisation of Peng et all 
(2010). The lines denote the best-fit model and its standard de- 
viation, according to Equation O 



(solid and dotted lines in Fig. [T2|. In principle, this signifi- 
cant anti-correlation allows us to use the colour-mass diagram 
as a diagnostic of the host properties when detailed observa- 
tional information, which would allow the determination of ac- 
curate star-formation rates and stellar masses, is not available, 
although the large scatter limits its reliability. Moreover, the 
AGN sample probed here is only a sub-sample of the total 
AGN population, selected both in the X-rays and at longer (far- 
infrared) wavelengths. There is evidence that the diflTerent selec- 
tions of AGNs bias their position in the colour-magnitude dia- 
gram (iHickox et al.Ll2009l) . with the X-ray- selected AGNs being 
in the green valley {{U - B) - (U - 5)rs ~ in Fig. [HI), the 
infrared- selected in the blue cloud and the radio- selected in the 
red sequence. The complex selection of the AGNs in our sample 
limits its representativeness. 



5. Discussion 

5.1. Is there an AGN-host correlation? 
5.1.1. Low redshifts (z < 1) 

Previous studies in the low-redshift Universe (iNetzeil 120091: 
ISerjeant & Hatziminaogloui l2009l) have shown a correlation be- 
tween the AGN luminosity (bolometric or optical) and the host 
galaxy luminosities (or star-formation rates), probing luminous 
QSOs with Lx > 10^^ erg s"^ and optically- selected QSOs, re- 
spectively. On the basis of our low-redshift (z < 1) - low- 
luminosity (Lx < 10"^^-^ erg s"^) AGN sample, we do not see 



sSFRms(z) 



(-1.36±0.17)[(t/-5)-(f/-5)Rs]+(0.2±0.4)(5) ^^^j^ ^i^^^ ^^^^^ ^ correlation between the AGN and the host 
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galaxy activity (parametrised by the sSFR; see Figs. [8] and \TU\ 
and Tables [T] and O. This implies that lower luminosity AGN 
activity, especially at low redshifts, is not directly linked to the 
state of the host galaxy, if the latter is parametrised by its sSFR 
or its "starburstiness". In the redshift range of the first three bins 
(z < 1.12) there are seven quiescent AGN hosts and 11 star- 
bursts, while the mean sSFR is within the borders of the main 
sequence. The selection of objects which have a FIR counter- 
part could aff'ect the mean "starburstiness" of our sample, but 
this eff'ect is found to be minimal when limiting the sample to 
the area where upper FIR limits are available. It is likely that 
the AGN process takes place as a result of instabilities that af- 
fect the nucleus but do not have any prominent eff'ect overall, 
being only confined to circumnuclear star formation. Indeed, 
there is a positive correlation of the AGN po wer with the nu- 
clear s tar formation in local Seyfert-1 g alaxies ([Thompson et al.L 
120091: iDiamond-Stanic & RieM 12012*), an d such a cor relation 
is also supported by recent models ( Hop kins & Quataeri.[201Q) . 
The causal mechanism behind this connection could be high- 
mass stellar winds fueling the AGN (see e.g. simulations by 
[Schartmann et al., 2009). This nuclear correlation however does 
not leave a clear mark on the overall observable properties of the 
system (infrared and X-ray luminosities). 



5.1.2. Higher redshifts {z > 1) 

In the more distant Universe there are studies finding both 
a correlation between t he AGN pow er and star-formation 
intensity dT richas et ap, 120091: iHatzi minaoglou et al. 201(1 
'Bonfield et al., 201 1), and no signs of any ( Seymour et al., 201 1; 
iRosario et al., .2012) . using diff'erent diagnostics and source se- 
lections, while there is evidence that the star-formation rates of 
the AGN ho sts are enhanced with respect to those at z < 1 
(iMullaney e t al., 2010). Combining data at diff'erent luminosities 
and redshifts. ILutz et al.l (2010) and lShao et al.1 (120 lOl) propose 
diff'erent mechanisms for the fuelling of both the star formation 
and the AGN, merger-driven for high luminosities and secular 
for lower, with the "high-luminosity" limit being strongly de- 
pendent on redshift ( see also Serjeant & Hatziminaoslou, 2009; 
IWilman et aH l2010h . Recently. iMullanev et al.1 (1201 2ah using 
deep Chandra and Herschel observations find that the increase 
of both infrared and X-ray luminosities with redshift aff'ect the 
observed Lx - SFR correlation, and do not detect it for moderate 
luminosity AGNs (Lx = 10^^-10^^ erg s"^) in individual redshift 
bins. However, such a correlation emerges at redshifts z ~ 1-2, 
if the stacked signal from individually X-ray undetected AGNs 
is factored in, rev ealing a similar Lx - relation to the SFR 
mam sequence (Mullaney et aU|2012b|). There is also evidence 
that the correlatio n is weaker at higher redshifts for the highest 
luminosity AGNs (iRosario et al.Ll2012h . 

Here, we use a combination of the deepest XMM-Newton 
and Herschel observations in combination with an SED de- 
composition techniq ue to remove any AGN flux from the far- 
IR wavelengths (see iMuUaney et al.L l201ll) . and find a corre- 
lation between the specific SF R and the X-ray lumin osity for 
z > 1 andLx > 10^^ erg s'^. 'MuUaney et al.' ('2012a") usmg a 
similar sample with a somewhat lower luminosity range (Lx = 
10^^ - 10^"^ erg s"^) fail to detect a significant correlation at those 
redshifts, suggesting that the higher luminosity sources are re- 
sponsible of the correlation detected in this work. Our sample is 
highly incomplete for Lx < 10^^ erg s"^ at z > 1 (see Fig. [6]). 
Following the discussion of the previous section, this behaviour 
can be either because the nuclear star formation in high X-ray 
luminosity objects is so strong that it dominates over that of the 



host, or there is a link between the AGN activity and the evo- 
lution of the host galaxy at higher redshifts and luminosities, 
parametrised by the overall sSFR. In the former case we would 
expect a correlation between t he obscuration of the A GN and 
the star-formation rate (see e.g.l Ballantyne et alll2006l) . at least 
for high luminosity objects showing some degree of obscura- 
tion, because of the expected increase in the covering factor and 
the column density of the obscuring material, if it is also re- 
sponsible for the star formation. In Figs. lllbl and lllclwe do not 
detect any such correlation; moreover, according to Ballantynd 
( 2008), a circumnuclear star-forming disk could not sustain very 
high star-formation rates (> 10MQyr"\ typical of rates of our 
sample) and is not compatible with a high-luminosity AGN, be- 
cause it would limit the necessary gas supply. This is an indi- 
cation that the star-formation rate is not nuclear and therefore 
not directly connected to the AGN obscuration. The lack of any 
correlation between the AGN obscuration and the sSFR also in- 
dicates that the star-forming gas is not directly connected to the 
AGN obscuration, so the obscu ration from the host galaxy (see 
iMartinez-S ansigr e et al.L 12009^ is not dominant. The apparent 
connection between the galaxy and AGN activity is therefore 
likely evolutionary. 

This correlation of star formation at galaxy scales with the 
AGN activity seems to be in disagreement with models sug- 
gesting that the AGN outflows quench the star-forming activ- 
ity by disrupting the cold gas supply (Di Matteo et al., 2008), as 
there are a number of AGNs with high luminosities, L2-iokeV > 
lO'^'^ergs"^ which are actively star-forming, with sSFR > 

1 Gyr"\ and the most active AGNs appear to be more "star- 
bursty" than lower L2-iokeV sources, at least in the redshift range 

2 < z < 3 (see Figs. [6] and [9b]). Such behaviour is consis- 
tent with the suggestion that the AGN activity might enhance 
the star-formin g activity o f the h ost galaxy instead of quench- 
ing it (see e.g. lElbaz et all l2009l) and one means of doing that 
is through the disruption of the density profi l e of t he host by 
an AGN-generated jet (see also iGaibler et al.L l201lh . An issue 
that has to be addressed in this case is whether the jet would 
be detected at radio wavelengths, since only 2/7 of the high- 
est sSFR and highest L2-iokeV sources in Fig. [8] are detected in 
the radio, one of them being only marginally radio-loud with 
^1.4 GHz = 6 X 10^^ erg s"^ Hz"^ We note that the radio luminos- 
ity of HE 0450 - 2958 (the source studied in lElbaz et al.ll2009h 
is not radio-loud according to our classification, having a radio 
luminosity in the 1.4-GHz band in the order of 10^^ erg s"^ Hz"^ 
(Feain et al., 2007), meaning that a relatively low radio luminos- 
ity jet could cause a star-formation episode. 

5.2. Where do AGNs live? 



Recent studie s (iDaddi et al.l. I2007L 120091: IPunne et al.1 
iPanneUa et all 120091: iMagdis et all l2oT( 



Elbaz et al.l. 



2009t 

2om 



have found a relation between the star-formation rate and the 
stellar mass, consistent with being linear at all redshifts from 
local to z ~ 4, but where the normalisation of this relation is 
strongly dependent on redshift (iKarim et al.L 120111: lElbaz et all 
2011). In this discussio n we use the star-formation "main- 
sequence" of Elbaz et al] (|201l h up to z = 2.156 and a constant 
value of sSFRms = 2 thereafter. As we can see in Fig. [9al the 
sSFRs of the AGN hosts are mostly on the main sequence, indi- 
cated by the grey area or above it. This is more clearly demon- 
strated in Fig. [9bJ where we plot the deviation from the main 
sequence ("starburstiness") of the AGN hosts. The black data- 
points and line denote the running mean (and the respective er- 
ror) of the whole AGN sample described in § 13.31 the line is 
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constant with redshift (within the errors) and close to the up- 
per border of the main sequence. This result also holds if we use 
a sample unbiased by the lack of upper limits for all the FIR- 
undetected sources (dashed line - see § 14.21) . This is a similar re- 
sult to XueetaL C201Q) who find that the SFR of AGN hosts is 
similar to that of non-AGN galaxies when using mass-matched 
samples for z < 3. 

Overall, there are 11 quiescent, 54 starburst and 34 main- 
sequence AGN host s, which is in agreement with the findings of 
ISantini et~aD (l2Q12h who use similar methods on a wider sam- 
ple. Within the luminosity range 10"^^ < Lx < lO'^'^ergs"^ we 
find 23/69 main-sequence, 38/69 starburst and 8/69 quiescent 
hosts (assuming an 1 cr confidence interval of a binomial dis- 
tribution). These numbers do not a gree at first glance with the 
findings of iMullaney et al.l (l2012ah w ho use a sample similar 
to the one used in this study. However, iMuUaney et al.l (l2012ah 
use a wider main- sequence region (a factor of three instead of a 
factor of two of the main-sequence sSFR) and if we adopt this 
definition, the above numb ers become 39/69, 27/6 9, and 3/69, 
respectively, much closer to iMuUaney et al.l (l2Q12ah . The resid- 
ual diff'erence of the fewer quies cent hosts found here is because 
of the stacking analysis done in iMuUaney et al.l (1201 2ah to esti- 
mate the behaviour of FIR undetected AGNs. The limited num- 
ber of sources in the "complete" sample in the area covered by 
Herschel-FACS d oes not allow us to perform such an analysis 
here. We note that lSantini et al.l (l2012h find similar results when 
they factor-in their stacking analysis of undetected AGNs. The 
increased mean sSFR of the AGN hosts we find in this study is 
in line with the Lx - sSFR correlation, suggesting that the AGN 
and star-formation processes are connected, either aff'ecting each 
other, or having a common cause. The most luminous AGNs 
with Lx > 10"^"^ erg s"^ (X-ray QSOs) are represented with green 
symbols in Figs.[9a|and[9bl and reside on average in the region of 
starburst galaxies (defined as having sSFR/sSFRms(^) > 2) for 
z > 2, which reflects the overall correlation between the AGN 
luminosity and the host activity. 

In the redshift range 1 < z < 2 there are a few high- 
luminosity AGNs which have very low sSFR and "starbursti- 
ness" values, placing them in the main sequence or even in 
the quiescent region. Although these objects are not enough to 
disrupt the sSFR-Lx correlation at those redshifts (see Figure 
[8]), they could be examples of the powerful AGN suppress- 
ing the star-formation. In a recent study, using Chandra X-ray 
data and Hersche l-SFIRE sub-mm (250 yum) data in the CDF-N, 
iPage et afl (l20lA find that the highest X-ray luminosity (Lx > 
lO'^'^ergs"^) AGNs are rarely detected in the sub- mm wave- 
length s, and therefore have modest SFRs (see also lTrichas et al.L 
l2012 h. In our sample, most of the X-ray QSOs (14/20 of the 
"complete" sample) are detected in the far-infrared, although at 
a shor ter wavelength (lOOyum) than in the sample of 'Page et al.' 
(l2012h . This could imply that there is some residual contribu- 
tion from the AGN in shorter FIR wavelengths. However, with 
our SED analysis we identify and remove the contribution of 
AGN flux i n the far-infra red flux, so this explanation is un- 
likely. Rosa rio et al.l (|201^ argue that the SFR-Lx relation starts 
to weaken above z ^ 1, and indeed the correlation we find is 
not very strong for the 1.156 < z < 1.599 redshift bin, as a re- 
sult of the low-sSFR QSOs in that redshift bin. We do find on 
the other hand, that at higher redshifts the sSFR-Lx correlation 
is stronger, and the high X-ray luminosity AGNs are on aver- 
age more "starbursty" than the overall sample. This could be a 
result of higher abundance of molecular gas at hig her redshift 
(see e.g. iDaddi et al.L 120101: iBournaud etall l201ll) . where de- 



spite the feedback from the powerful AGN, the star-formation is 
still powerful. 

6. Conclusions 

We select 131 AGNs from the 3 Ms XMM-Newton survey and 
measure their star-formation rates using long wavelength far-IR 
and sub-mm fluxes with rest- frame wavelength above 20 jim. For 
32 of the 131 sources we are able to derive only an upper limit 
of the star-formation rate. We take special care in modelling the 
spectral energy distributions, identifying and removing the AGN 
contribution, and derive the sSFR and stellar masses of the hosts, 
comparing them to the AGN properties (X-ray luminosity and 
absorption). Our results can be summarised as follows: 

1 . We find no evidence for a correlation between the sSFR and 
the X-ray luminosity for sources with Lx < 10^^-^ erg s"^ and 
atz < 1. 

2. We find a correlation between the sSFR and the X-ray lumi- 
nosity for sources with Lx > 10^^ erg s"^ and at z > 1. There 
is no indication that this correlation is a result of a redshift 
eff'ect, as it is present even when we divide the data into nar- 
row redshift bins. We argue that it is instead a result of the 
AGN-host co-evolution, which is more prominent for higher 
luminosity systems, confirming previous results. 

3. We do not find any correlation between the star-formation 
rate (or the specific SFR, or the "starburstiness") and the 
X-ray absorption derived from high-quality XMM-Newton 
spectra, at any redshift or X-ray luminosity. We assume that 
this is an indication that the X-ray absorption is linked to the 
nuclear region, and the star-formation to the host. 

4. Comparing the sSFR of the hosts to the characteristic sSFR 
of star-forming galaxies at the same redshift ("main se- 
quence") we find that the AGNs reside mostly in main- 
sequence and starburst galaxies, with the mean specific SFR 
being close the limit between main-sequence and starburst 
hosts. This reflects the AGN-starburst connection. 

5. Higher X-ray luminosity AGNs (X-ray QSOs with Lx > 
10"^"^ erg s"^) are found in starburst hosts with average sSFR 
more than double that of the "main sequence" at any redshift 
above z ~ 2. At lower redshifts (z ~ 1.5) we find a num- 
ber of QSOs with low sSFR values, which drive the mean 
starburstiness of QSOs to a value consistent with that of the 
overall AGN population. 

6. We test the reliability of the colour-magnitude diagram in 
assessing the host properties, and find a significant anti- 
correlation between the "redness" (deviation of the rest- 
frame colours from the line dividing red and blue galax- 
ies, without any correction for AGN contribution or dust ex- 
tinction), and the "starburstiness" (the sSFR divided by the 
"main sequence" sSFR at a given redshift). 
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